A CONTINUOUS MODEL OF TRANSPORTATION REVISITED 



LORENZO BRASCO AND MIRCEA PETRACHE 



Abstract. We review two models of optimal transport, where congestion effects during 
the transport can be possibly taken into account. The first model is Beckmann's one, 
where the transport activities are modeled by vector fields with given divergence. The 
second one is the model by Carlier et al. (SIAM J Control Optim 47: 1330-1350, 2008), 
which in turn is the continuous reformulation of Wardrop's model on graphs. We discuss 
the extensions of these models to their natural functional analytic setting and show that 
they are indeed equivalent, using an ad hoc generalization of Smirnov decomposition 
theorem for flat 1— currents, which may fail to be normal. 
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1. Introduction 

1.1. Theoretical background. The present work is motivated by the study of transport 
problems for distributions, started in [6] and [12], which we want to try and connect to 
related works in the theory of currents present in [TJ [27] [36] . A motivation for such an 
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extension is the study of distributions of the form 

oo 

^(<*Pi - <*Qi)> with ^2\Pi - Qi\ < oo, 

i=l i=l 

as in [33]. Such distributions arise as topological singularities in several geometric varia- 
tional problems as described for example in [HJ [23j [32J, l34l 135]. 

To start with, we formally define three variational problems which can be settled (for 
simplicity) on the closure of an open convex subset Q C R , having smooth boundary. 
For the moment, we will be a little bit imprecise about the datum /, but we will properly 
settle our hypotheses later. The first problem is the minimization of the total variation of 
a Radon vector measure, under a divergence constraint: 



(B) min <^ / d\V\ : -divV = f, V ■ u n = 

v [Jn 

The above problem can be connected by duality with the following one, called the Kan- 
torovich problem: 

(£) max {(/,</>> : ||V0[| L oo (n ) < l} , 

where now the variable (ft is a Lipschitz function and (•, •) represents a suitable duality 
pairing. Finally, the third problem is the minimization of the total length 

(M) mm|y ^( 7 )dQ( 7 ) : ( eo - ei ) # Q = /|, 

where V is the space of Lipschitz continuous paths 7 : [0,1] — > Q, the length functional I is 
defined by 



eo, e\ are the evaluation functions giving the starting and ending points of a path, and the 
variable Q is a measure concentrated on V . 

The classical setting for the above problems is when / is of the form / = f + — f~ , where 
/ + and /~ are positive measures on having the same mass (for example, conventionally 
one can consider them to be probability measures). We point out that in this case, a more 
familiar formulation of is certainly the so-called Monge-Kantorovich problem 



(M') mini / \x - y\ drj(x, y) : (tt x ) # t] = f + and (%)#?? = f 
*i Unxn 

where ir x ,TTy : 17 x 0, — > stand for the projections on the first and second variable, 
respectively. It is useful to recall that the link between (I.Mj) and &A4'\i is given by the fact 
that if 770 is optimal for Monge-Kantorovich problem then the measure which concentrates 
on transport rays, i.e. 

Qo = Sxy drjo(x , y) , where ~xy stands for the segment connecting x and y, 
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is optimal in ()7WP and 



Kl) d Qo(l)= \x-y\drj (x,y). 
v Jnxn 



When / has the above mentioned form f + — f~, the equivalence of the three problems above 
is well understood: the equivalence of (|7W / |) = (|A^P and (JQ is the classical Kantorovich 
duality (see [24]), while that between (JHJ) and (IXH) seems to have been first identified in 

m- ' 

Recently the equivalence of the above three problems has been shown in [6j for / belong- 
ing to a wider class, i.e. when / is in the completion of the space of zero-average measures 
with respect to the norm dual to the C 1 (or flat) norm. This wider space was studied in 
[22j and characterized recently in [HIT]. A different point of view is also available in [25J, 
where the space of such / is called W > . 

1.2. Goals of the paper. Our starting observation is that problem pertains to a wide 
class of optimal transport problem, introduced by Martin J. Beckmann in [2], which are of 
the form 

(Bu) min | J H{V) dx : divV = f, V ■ vq = (A , where c x \z\ p < H(z) < c 2 \z\ p , 

for a suitable density-cost convex function T~L : — > M + and p > 1. For a problem of 
this type, the question of finding equivalent formulations of the form (jAp) and (|A4p has 
already been addressed in [10] (see also [8]), under some restrictive assumptions on /, like 
for example 

/ = /+-/- with /+, /- € LP(n) and //+=//" = 0. 

Jn Jn 

The goal of this paper is to complement and refine this analysis, first of all by studying 



problem <\Bfi\) in its natural functional analytic setting, i.e. when / belongs to some dual 
Sobolev space W~ l ' p (whose elements are not measures, in general). Also, by expanding 
the analysis in [H [10] , we will see that alternative formulations of the type ([AC]) and (|A4p 
are still possible for (\Bu\) in this extended setting. These formulations are still well-posed 
on the dual space W~ l ' p and equivalence can be proved in this larger space. The problem 
corresponding to (JQ will now have the form (see Section 3 for more details) 



(/C%) max 



(/,</>>- / n*(v^)dx 



n 



and the equivalence with ( By ) will just follow by standard convex duality arguments (that 
we repeat in this paper, for reader's convenience). On the contrary, in the proof of the 
equivalence between ( By ) and its Lagrangian formulation 



(M n ) nun |_£^(^q) dx : ( e o - ei)#<9 = / J 
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some care is needed, due to the fact that / now is not a measure. Here the measure iq will 
be some sort of transport densit$\ generated by Q, which takes into account the amount 
of work generated in each region by our distribution of curves Q (see Section U] for the 
precise definition). In particular, the proof of this equivalence will point out another not 
emphasized connection to Geometric Measure Theory. 

The main result of this paper can be formulated as follows (see Theorems 13.11 and 14.51 
for more precise statements): 

Main Theorem. Let 1 < p < oo. Suppose Q C is the closure of a smooth bounded 
open set, let f E W~ l,p (VL) and let % be a strictly convex function having p— growth, with 
p > 1. Then the minima in (By), (My) an d the maximum in (JCy) are achieved and 
coincide. Moreover we have the following relationship among the optimizers of the three 
problems: 

(i) the unique minimizer of (By) corresponds to a minimizer of (My.), in the sense 
of Proposition 

(ii) each minimizer of (My) corresponds to the unique minimizer of (By), in the sense 
of Proposition 

(iii) the unique minimizer V of (By) and the unique maximizer v of (By) are linked by 
the relation V{x) 6 d~H*(x,Vv(x)) a.e. in Vt, as specified in Theorem \3.1{ 



The connection of the first two points in the above theorem to Geometric Measure 
Theory lies in the basic theory of flat 1- currents, of which we recall the first steps in the 
(long) appendix at the end of the paper. Indeed, in order to show equivalence of (\By\j and 
A-M'hD i n full generality, i.e. in the space of distributions W~ 1,p , the cornerstone will be 
Smirnov decomposition theorem for 1— currents. We also notice that to this aim we need 
to slightly extend Smirnov's result to cover the case of flat currents with finite mass, i.e. 
we drop the assumption on the finiteness of the mass of their boundaries (see Theorem 
IA,20[) . Translating this in the language of vector fields V, this extension is needed since 
the divergence of V now will be just an element of W~ 1,p , rather than a measure. 

For the sake of completeness and in order to neatly motivate the studies performed in 
this paper, it is worth recalling that the proof of this equivalence in |10j was based on 
the Dacorogna-Moser construction to produce transport maps (see [13]), which has been 
revealed a powerful tool for optimal transport problems^- In a nutshell, this method consists 



in associating to the "static" vector field V admissible in < \By\ ), the following dynamical 
system 

din = div ( {1 _ t) J + + tf _ tn) , Mo = / + , 



^When H(t) — \t\, problem JA^hI I is again the Monge-Kantorovich one and iq for an optimal Q is 
nothing but the usual concept of transport density, see 1151 120j . 

It is worth remarking that the first proof of the existence of an optimal transport map for problem 
(f^7J> more than 200 years after Monge stated it, was based on a clever use of this construction (see |18|). 
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i.e. a continuity equation with driving velocity field Vt given by V rescaled by the linear 
interpolation between f + and /~. Supposing that one can give a sense (either deterministic 
or probabilistic) to the flow of Vt, then the construction of the measure Qy concentrated 
on the flow lines of Vt paves the way to the equivalence between the Lagrangian model 
flA^D and flff-^D (see [U [10] for more details). However, this strategy seems to need some 



restriction on the datum /, in particular / should be a Radon measure on Q. 

1.3. Plan of the paper. In Section[2]we describe the function spaces W~ 1,P (Q) and prove 
the existence of a minimizer for (Bu)- In Section [3] we prove the equivalence of {By) with 
{JOh), by appealing to classical convex analysis formulas. The aim Section[4]is to introduce 
the Lagrangian counterpart of Beckmann model and to show how the two models turn out 
to be equivalent. A self-contained Appendix complements the paper. There we introduce 
relevant concepts from Geometric Measure Theory and prove the extension of Smirnov's 
decomposition theorem to general flat 1— currents with finite mass (Theorem IA.20j) . which 
is the main ingredient for the proofs of Section [H 

2. Well-posedness of Beckmann problem 

Let Q C M N be the closure of an open bounded connected set, having smooth boundary. 
In particular, in what follows f2 will always be compact. Given 1 < q < oo, we indicate 
with W^ q {9) the usual Sobolev space of L q {Q) functions, whose distributional gradient is 
in L q {Q) as well. We then define the quotient space 

where ~ is the equivalence relation defined by 

u ~ v there exists c G R such that u(x) — v{x) = c € M. for a.e. 

When needed, the elements of W 1,q {£l) will be identified with functions in W l,q {£l) having 
zero mean. We endow the space W 1,q (yi) with the norm 



ulU,.„^ ■=( f Q \vu(x)\ q dx^j 9 , tie w^ q (n), 



then we denote by W 1,p (f2) its dual space, equipped with the dual norm, defined as usual 
by 

\\ T \\w-i,r(Q) ■■= SUp{(T,(/>) : <p G W 1 *^), \\<p\\ Wl , q = l} , 
where p = qj '{q — 1). We start recalling the following basic fact. 
Lemma 2.1. Let T G W~ l > p (n), then 



max \(T,(p)\-~ [ \Vip\ q alx 
<peW 1 'i(Q) q Jq 
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Proof. For every 93 G W 1,q (£l), we clearly have 

\{T,tp)\-- [ |V^| 9 da;<sup A|(I»|- — I \Vip\ q dx 
Q Jn a>o L Q Jn 

On the other hand, the supremum on the right is readily computed: this corresponds to 
the choice 



\V(f\ q dx 



1 

"9-1 



which gives 



sup 

A>0 



A|(7»|- — / \Vu\ q dx 
Q Jn 



1 f\(T,<p)\Y 



Passing to the supremum over tp G W 1,q (Q,) and using the definition of the dual norm, we 
get the thesis. □ 

We also denote by £[(Q) the space of distributions of order 1, with (compact) support 
in Q. In what follows, we will tacitly identify this space with the dual space of the Banach 
space C 1 (r2), endowed with the norm 

IMIc 1 ^) = IMU°°(n) + ||V<^||ioc»(n)- 

By denoting i/q the outer normal versor to <9f2, we have the following characterization for 
the dual space W~ l ' p {Q). 

Lemma 2.2. Let p = q' = q/{q - 1). Given a vector field V G D>(Sl;M. N ) and T G £[{VL), 
we say that V satisfies 



(2.1) 

if 



— div V = T in $7, V ■ vq = on <9$7, 
V(f(x) ■ V(x) dx = (T, tp), for every tp G C 1 (Q). 



// we set 

£l p (Q) = {T G £[(n) : there exists V G L p (n-,R N ) satisfying (gTQ }, 
we then have the identification 

w- 1 *(n) = £l P (n). 

Proof. Let T G W~ 1,P (Q), first of all we observe that then T G £\{Q) as well. Now, consider 
the following maximization problem 



sup (T,v) - - J \Vv{ 
«eH"-«(!i) 1 



x)\ q dx. 
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By means of the Direct Methods, it is not difficult to see that there exists a (unique) 
maximizer u € W l ' q (£l) for this problem. Moreover, such a maximizer satisfies the relevant 
Euler-Lagrange equation, given by 

/ \Vu{x)\ q ~ 2 Vu(x) ■ Vip(x)dx = {T,(p), for every ip G W hq (n). 

By defining V = \Vu\ q ~~ 2 Vu € L P (H;M. N ), the previous identity implies T € f( p (0). 
Conversely, let us take T € £{ (fi), then for every tp G C 1 (r2) equation (|2.1|) implies 



l<7»| 



Vip(x) ■ V(x) dx 



< 



\\ V \\lp(Q)- 



Using the density of C 1 (f]) in W '^(O,), we then get that T can be extended in a unique 
way as an element (that we still denote T, for simplicity) of W ' q (Q,). Observe that this 
extension satisfies 

ll T llw- 1 .i'(n) - H^lk^n)) 
by taking the supremum in the previous inequality. □ 

Remark 2.3. We remark that the elements of £[ p (Q) have "zero average", i.e. 

(T, 1} = 0, 



as follows by testing the weak formulation of (12. ip with <p = 1. This is coherent with the 
previous identification W~ l,p {£l) = £[ (Q), since by construction the space W 1,q {£l) does 
not contain any non trivial constant function. 

Example 2.4. Consider the measure T = 5 a — 5b for two points a ^ b € M. N . We claim 
that 

T = 5 a -5 b eW~ 1 ' p (n) if and only if 1 < p < N/(N - 1), 

where O is a sufficiently large ball containing a, b in its interior. We will prove this by using 
the characterization of Lemma [23 



Suppose indeed that there exists some V E L p (il), such that — divU = T. We pick a ball 
B r (a) centered at a and having radius r such that 2r < \a — b\. Then given e < r, we take 
a Cq{B t {o)) function r\ e such that 

r} e = \ in £? r _ e (a) and ||V%||l°° < Ce" 1 . 
Thanks to our assumption, we will have 

1 = (T,rj E ) = f V ■ Vrj E dx, 

JB r {a) 

so that 



I \V\dx>^. 
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1 

r 



By Holder inequality, this easily implies a lower bound on the L p norm of V, namely 

i-p 

\V\ p dx > e p \B r {a) \ B r _ E (a)\ l -P = C N)P e p r N ^^ 1 

'B r {a) 

We now make the choice e = r/2, so that from the previous we can infer 

\V\ p dx > C N;P r p+N ^ = Cn^-p^-V. 

B r (a) 

The previous estimate clearly contradicts the assumption V £ L p (£l), if the exponent 
N — p (N — 1) is not strictly positive. Therefore we see by Lemma 12.21 that p < N/(N — 1) 
is a necessary condition for T € W' 1 ^^!). 

This condition on p is also sufficient for T to belong to as we will now show. Let 

us set 2 r = \a — b\, for simplicity let us suppose that a = (— r, 0, . . . , 0) and b = (r, 0, . . . , 0). 

we 



Then, using the notation x = (xi,x') for a generic point in M. , where x' E 
consider the following vector field 



V a ,b{x) 



(Xl+T, X 1 ) 

(x 1 + r) N 

(Xl - T,X') 



[XI - T 



t (o,...,o), 



if Ix'l < r and \x'\ — r < x\ < 0, 



if \x'\ < t and r — \x'\ > x\ > 0, 



otherwise. 



It is easily seen that div V a ,b = S a — 5b and that V a ,b is supported on the set 



D a ,b = < {xi,x') £ 



> \x I + \xi 



which is just the the union of two cones centered at a and b, having opening 1 and height 
t = \a — b\/2. Also, by construction we have 



\V a ,b{x)\ P dx 



Qa 



+ \x 



l\2 



2 E ¥Nu N I (X! + T )-Np+P+N-1 ^ 



t J {x' : \x'\=xi+t} 
■0 



(xi + r) Np 



dx dx\ 



so that finally 

ii^iii^c^ia-^r-^- 1 ), 

thanks to our assumption p < N/(N — 1). For some related constructions, the reader is 
referred to [H Proposition 3.2] and [32l Lemma 8.3]. 

As a consequence of Lemma 12.21 we have the following well-posedness result for Beck- 
mann problem. 
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Proposition 2.5. Let H : ft x M. N be a Caratheodory function, such that z i— > T~L{x,z) is 
convex on M N , for every x G fl. We further suppose that % satisfies the growth conditions 

(2.2) \{\z\ p - 1) < H(x,z) < \(\z\ p + l), (x,z)£nxR N , 

A 

for some < A < 1. Then the following problem 

(2.3) min j / U(x,V)dx : -divV = T, V ■ u n = 0} , 

V€LP(Q;R N ) Utl " J 

admits a minimizer with finite energy if and only ifT G Ty _1,p (f2). 

Proof. Let T G H /_1,p (0), then thanks to Lemma 12.21 there exists at least one admissible 
vector field Vq with finite energy, so that the infimum (|2.3p is finite. If {V^} ne N C L p (fl; M N ) 
is a minimizing sequence, then the hypothesis on H guarantees that this sequence is weakly 
convergent to some V G L P (Q; M. N ). Thanks to the convexity of T~L, the functional is weakly 
lower semicontinuous, i.e. 

H(x,V)dx < liminf / H(x,V n )dx 
n ^°° Jn 

= min < / T-L(x,V(x))dx : T f* V ^ \- 
VeLP(Q;R N ) [Jq V-vq = Q J 

Moreover, this limit vector field V is still admissible, since 

/ V<f ■ V dx = lim / V</9 ■ V n dx = (T, ip), for every <p G C 1 (fl), 
Jn n ^°° Jn 

by weak convergence. Therefore V realizes the minimum. 

On the other hand, suppose that T G" W~ 1,p (Vl). Again thanks to Lemma [2. 2 1 we have that 
the set of admissible vector fields is empty, so the problem is not well-posed. □ 

For the sequel, we need the following definition. 

Definition 2.6. We say that a vector field V G L} (fi;R N ) is acyclic if, whenever we can 
write V = V\ +V%, with \V\ = \V1\-\-\V2\ and div V\ = in distributional sense, there must 
result Vi = 0. 



The following is a mild regularity result for optimizers of (12.3f) in the isotropic case, i.e. 
when % depends on the variable z only through its modulus. This will be crucial in order 
to equivalently reformulate (|2.3|) as a Lagrangian problem, where the transport is described 
by measures on paths. 

Proposition 2.7. Let us suppose that Tl satisfies the hypotheses of Proposition \2.5[ In 
addition, we suppose that 

z 1 — y ?i(x, z) is a strictly convex increasing function of \z\, for every x. 

Then there exists a unique minimizer V for f|2.3|) and V is acyclic. 
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Proof. The uniqueness of V just follows by strict convexity, so let us prove that V is acylic. 
Let us suppose that we can write V = V\ + V2, for some vector fields V\, V2 G L 1 (f2;lR Ar ) 
such that 

\V\ = |Vi| + \V 2 \ and divyi = 0. 

As a consequence, we have divF = divV2 and \V\ > |T^|. Thus V2 is a competitor for 
problem (j2.3|) with energy not larger than that of V, thanks to the monotonicity of Tl. 
Since V is the unique minimizer, it must have energy equal to that of V2. Thus \V\ = | V2 1 
and \Vi \ =0 almost everywhere. Owing to the definition, this shows that V is acylic, thus 
concluding the proof. □ 

3. Duality for Beckmann problem 

In what follows, we will need the following general convex duality result (for the proof, 
the reader is referred to |17t Proposition 5, page 89]). The statement has been slightly 
simplified, in order to be directly adapted to our setting. 

Convex duality. Let T : Y — > R be a convex lower semicontinuous functional on the 
reflexive Banach space Y . Let X be another reflexive Banach space and A : X — > Y a 
bounded linear operator, with adjoint operator A* : Y* — > X* . Then we have 

(3.1) sup{x*,x}-J r (Ax)= mm{T*(y*) : A*y* = x*}, x* E X* , 

xeX y*&* 

where J 7 * : Y* ->1U {+00} denotes the Legendre-Fenchel transform of J 7 . Moreover, if 
the supremum in (|3.ip is attained at some xq G X , then the minimum in (|3.ip is attained 
as well, by a i/q G Y* such that 

y* G d-F(Axo). 

Thanks to the previous result, we obtain that Beckmann problem admits a dual formu- 
lation, which is a classical elliptic problem in Calculus of Variations. 



Theorem 3.1 (Duality). Let % be a function satisfying the hypotheses of Proposition W. 
andT G W~ 1 ' P (Q). Then 



mm 

VeLP(n-A N ) <- Jn 

(3.2) 



H(x, V(x)) dx : -div V = T, V ■ u n = o} 



max <(T,v)— / 7i*(x, Vv(x)) dx 



where H* is the partial Legendre-Fenchel transform ofTL, i.e. 

n*(x,0 = sup £-z-H(x,z), xeSl,£€R N . 



Moreover, if Vo G L p ($7) and vo G W 1,q {VL) are two optimizers for the problems in (|3,2p . 
we have the following primal-dual optimality condition 

(3.3) V £ dH*(x,Vv ) inQ, 
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where dTi* denotes the subgradient with respect to the £ variable, i.e. 

dH*(x,€) = {zeR N : H*(x,£) + U(x,z) = £ ■ z, for every x G 0}. 
Proof. To prove (|3.2p , it is sufficient to apply the previous result with the choices 

Y = L q (VL; R N ), X = W 1 ' q (n), F(<f>) = [ H*(x, <j>{x)) dx and A (tp) = Vtp. 

Jn 

Observe that the operator A is bounded, since 

P(y)l|y = l|V^|| L 9(n) = IMIx, for every tp G X, 

and that 

[ H**(x,£(x))dx = [ H(x,£(x))dx, 
Jn Jn 
since £ i— > ?i(x,£) is convex and lower semicontinuous, for every x G f2. We only need 
to compute the adjoint operator A* : L P (Q,;R N ) —> W~ 1,p (£l). Let us define the map 
D : L P (Q;R N ) — >• S[ p (U) by 

DV G S[(n) such that (DV, tp) = [ Vtp- V dx, for every <p G C x (fi). 

Jn 

Observe that this is a linear operator, whose image is contained in = W~ l *{n) 

just by construction and by definition of £i tP (£l). Moreover, for every tp G C 1 (r2) and 
V G L P (£1;R N ) we have 



(Ay>,y)= / Vtp(x) ■ V(x) dx = (tp, D V). 
Jn 

By density of C l (9) in W^ifl), we obtain that D = A*, then fl3jZ]) follows from (pTT) . 

The primal-dual optimality condition f|3.3|) as well is a direct consequence of the second 
part of the convex duality result. It is sufficient to observe that the maximum in (|3.2p 
is attained at some vq G W 1,p (Q), simply using the Direct Methods. This implies that a 
minimizer Vq of Beckmann problem has to satisfy 

V G dT(Vv ), 

which implies directly (|3,3p . □ 

A significant instance of the previous result corresponds to T~L(x,z) = \z\ p . Thanks to 
Lemma 12.14 we have the following. 

Corollary 3.2. For every T G W~ 1,p (d), we have 

W T \\w-i.p(n) = m ; n „, (ll^lli^) : -divV = T, V-uq=o\. 

Proof. It is sufficient to use (|3.2p and Lemma 12. 11 paying attention to the easy fact 

max \(T,tp)\-- I \Vtp\ q dx= max (T,<p) - - [ \Vtp\ q dx. 
V €W 1 -i(n) q Jn ^ew^^n) q Jn 

This establishes the thesis. □ 
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Corollary 3.3. Under the hypotheses of Theorem \3.1\ we have that the functional 

$ H : W-^{n) R + 

T h- > minimal value (12 .3p 

is convex and weakly lower semicontinuous. 

Proof. It is sufficient to observe that thanks to Theorem 13. 1\ the value (|2.3p can be written 
as a supremum of the linear functionals L v defined by 

L^T) = (T, ip) - [ H*(x, Vip) dx, <p G 
Jn 

Then the thesis follows. □ 

About the duality result of Theorem 13-H some comments are in order. 

Remark 3.4 (Economic interpretation). By the so-called Legendre reciprocity formula in 
Convex Analysis, the primal-dual optimality condition (|3.3|) can be equivalently written as 

(3.4) Vw € &H(x,V ), in ft, 

so this result is the rigorous justification of the necessary optimality conditions derived in 
[2j Lemma 2]. Such a function vq will be called a Beckmann potential and its economic 
interpretation is that of an efficiency price system. It can be seen as a generalization of 
a Kantorovich potential, to a situation where the cost to move some unit of mass from x 
to y is not fixed, but it depends on the quantity of traffic generated by the transport Vq 
itself. Heuristically, observe that in such a situation this minimal cost will be given by the 
"congested metric" 

d Vo (x ,xi)= min / | VH(-, Vq) o 7] \j'(t)\ dt. 

This means that each mass particle is charged for the marginal cost it produces, the latter 
being the derivative of the function % (we suppose for simplicity that H has a true gradient 
and not just a subgradient). Then vo acts as a Kantorovich potential for the Optimal 
Transport problem 



mm 



d Vo (x,y)dn(x,y) : (ir x ) # n = T + and (%)#77 = TL 



where we suppose for simplicity that T = T + — T_, with T + and T__ positive measures 
having the same mass. It should be remarked that V</?o does not give the direction of 
optimal transportation in Beckmann problem, since V^o an d Vo are linked through the 
relation ()3.4|) . thus they are not parallel in general. Actually, this is the case only when the 
cost function rl is isotropic, i.e. when it just depends on |V|, for every admissible vector 
field V. This is the case studied by Beckmann in his original paper [2]. 

Remark 3.5 (Regularity of optimal vector fields). We point out that if z 1— > H(x,z) is 
strictly convex, then £ i-> H*(x, £) is C . Thus in this case, the optimal Vq is unique and 
we have 

V Q = S/H*(x,Vv ). 
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Then the regularity of the optimal vector field Vq can be recovered from the regularity of 
a Beckmann potential, which solves the following elliptic boundary value problem 

-divVH*(x, Vu) = T, in fi 
V%*(x,Vu)-vq = 0, on da 

For instance, if H* in uniformly convex "at infinity" , i.e if there exists M > such that 



(3.5) 



1 + \z\ 2 )^ < Cmm(D 2 H*(x,z)£,£), for every \z\ >M,xeU, 
\D 2 H*(x,z)\ < (l + |z| 2 )^, (x,z) GflxK* 



and 



then V is bounded, provided that T G L N+£ (Cl), since in this case solutions to (13.51) are 
Lipschitz. These assumptions are verified for example by (see |10j ) 

H*(z) = \\z\-5) q + , zeR N , 

where (■)+ stands for the positive part and 5 > 0, but are violated by anisotropic functions 
of the type 

N 



n*(z) = J2-(\z i \-S i )l, ze 



1=1 * 

considered for example in [BJ J9| . 

4. A LAGRANGIAN REFORMULATION 

The aim of this section is to introduce a Lagrangian counterpart of Beckmann model and 
to show how the two models turn out to be equivalent. The model we are going to present is 
a continuous version of the classical discrete model by Wardrop, which we mentioned above. 
This continuous model has already been introduced in [12] and the equivalence has been 
discussed in [ID]. Here, paralleling the case of Beckmann model, we go one step further, by 
extending this model to its natural setting, i.e. the dual space W~ 1,P (Q). Moreover, we will 
prove equivalence of the models without imposing any additional regularity on the datum 
T. This will be achieved by means of Smirnov decomposition theorem for 1— currents (see 



Given two Lipschitz curves 71, 72 : [0, 1] — > fi, we say that they are equivalent if there exists 
a continuous surjective nondecreasing function t : [0, 1] — > [0, 1] such that 

72 (*) =7i(t(*)), for every t e [0, 1]. 

Let us call £(f2) the set of all equivalence classes of Lipschitz paths in U. We can introduce 
a topology on this set by defining the following distance 

^(71,72) : = rnax{|7i(i) - *jb(t)\ : t G [0, 1], ji equivalent to 7J . 

Observe that convergence in this metric is nothing but the usual uniform convergence, up 
to reparametrizations. 
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We denote the class of cr-finite positive Borel (with respect to the above topology) measures 
on £(Q) by ~#+(£(Q)). For Q G ^# + (£(f2)), we may define the corresponding traffic 
intensity by 

(i Q ,<p):= I ( [ tp(j(t))\ 7 '(t)\dt) dQ( 7 ), ^C(O), 

provided that the outer integrals converges, in which case we say that "the traffic intensity 
iq exists" . If this is the case, then the following integral 



(i Q , <p) = f ( f <p(r,(t)) ■ j'(t) dt) dQ( 7 ), <p G C(0; ] 
Jan) \Jo J 



also converges. These definitions do not depend on the particular representative of the 
equivalence class we chosen, since the integrals in brackets are invariant under time repa- 
rameterization. 

Remark 4.1. Observe that iq counts in a scalar way the traffic generated by Q, while 
ig computes it in a vectorial way. This means that in principle iq and \iq\ could be very 
different, since in iQ two huge amounts of mass going in opposite direction give rise to a lot 
of cancellations, as the orientation of curves is taken into account. As a simple example, 
suppose to have two distinct points xq ^ xi, then we consider the measure 

1 1 

Q — 2 7i "i" 2 72 ' 

with 7i(t) = (1 — t) xo + 1 x\ and o" 2 (t) = (1 — t) x\ + txo, that is we simply exchange the 
mass in xo with that in x\ and vice versa. By computing the traffic intensity, we get 

which takes into account the intuitive fact that on the segment xo^T globally there is a non 
negligible amount of transiting mass. On the other hand, it is easily seen that 

Given T G we define also the following space of measures 

Qp(T) := {Q G Jt+(C(Q)) : iq G L p (tt) and (e - ei) # Q = r}, 

where ej : £(17) — > £1 is defined by e^) = for i = 0, 1 and the equality (eo — &i)#Q = T 
has to be understood in the distributional sense, i.e. 

/ [99(7(0)) - p(7(l))] dQ(7) = <T, 95), for every ip G C 1 ^)- 
J£(Q) 

The convergence of the above integral is again in the sense of distributions. 

Remark 4.2. All these definitions are reformulated in Appendix |A] in terms of currents. 



A CONTINUOUS MODEL OF TRANSPORTATION REVISITED 15 

Let us now consider T~L : f2 x M + — > R + a Caratheodory function such that 

(4.1) X(t p -1) <H(x,t) < \(t p + l), xeQ,t€R + 

A 

for some < A < 1 and such that 

t i-> ~H(x,t) is convex, for every x£(]. 
If Q P (T) ^ 0, then we may define the following minimization problem: 



(M H ) inf / H{x,i Q (x))dx. 

QeQ P (T) J n 

We will show that the above minimization is well-posed and equivalent to the one in (|2,3p . 
At this aim, we will use the following result, proved in the appendix (see Theorem IA.20j) 
in a slightly stronger formulation in terms of currents. 

Proposition 4.3. Let 1 < p < oo. Suppose that V G L P (Q,M. N ) and that V is acyclic. 
If we set T = — div V , then it is possible to find a positive a— finite Borel measure Q € 
^ + (C(£l)) such that 

(eo — ei)#Q = T in distributional sense. 

Moreover, we have 

\q = V and iQ = \V\. 

In particular Q G Q p (T). 

Thanks to the previous result, we can at first give a necessary and sufficient condition 
for the set Q P (T) to be not empty. 

Proposition 4.4. The set Q P (T) is not empty if and only if T € W~ 1,P (Q). 

Proof. Let us suppose that T g" W~ 1,P (Q) and assume by contradiction that there exists 
Qo € QpiT). In particular 



(4.2) / \i Qo {x)\ p dx < +oo. 

Jn 

The vector measure iq satisfies (|2.ip . since 



f V V ■ di Qo = [ ( [ 
Jn Jc(n) \Jo 



V^(t))- 7 '(t)dtj dQ ( 7 ) 
p( 7 (l)) - p( 7 (0))l dQ ( 7 ) = (T^}, 

<c(n) ' 

for every ip € C 1 (J7). Also, thanks to the fact that |ig | < iQ and to (|4.2p . we have that 
i Qo G L p (Q;R N ). This contradicts the fact that T G" I^ _1 ' p (0), as wanted. 

Let us now take T G W~ 1,P (Q), then there exists a minimizer V of (12. 3D . Thanks to 
Proposition 12.71 we know that V is acyclic, so that by Proposition 14.31 we can infer the 
existence of Q G Q P (T). This gives directly the thesis. □ 
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We can now prove our equivalence result, which is the main result of this section. Observe 
that we prove at the same time existence of a minimizer for fl-M-^p . 

Theorem 4.5. Let Ti : Q x M + — > R + be a Caratheodory function satisfying (|4.ip and 
such that 

t h- > Ti(x,t) is convex and increasing, x € fi. 
For every T £ VF~ 1,p (il), we /iai>e 

(43) QefeV) 1 H ^ %Q{X)) dx = yjfcw {L % ^ ^ dX ~V Y "n = ' } ' 

and the infimum on the left-hand side is achieved. 

Moreover, if Qo £ Qp(T) is optimal, then iq € L P (£1;~R N ) is a minimizer of Beckmann 
problem. Conversely, ifVo is optimal, then there exists Qy £ Q P (T) such thatiQ Vo = \iq Vo \ 
minimizes the Lagrangian problem. 

Proof. Since T £ W~ 1,p {£i), we first observe that by the previous result the set Q P (T) is 
not empty. Moreover, for every admissible Q, we have \'iq\ < iq and \q is admissible for 
Beckmann problem. Using the monotonicity of ~H(x, •), we then obtain 

min < / H(x,\V\) dx : T f* V ^ !F' \ < inf / %{x, in(x)) dx < +oo. 
veLp(n;RN)\Ja V-vq = J Q&Q p (T)J n 

Let us now take a minimizer Vo £ L^f^R^) for Beckmann problem, then by Proposition 
12.71 this is acyclic. Thus by Proposition 14. 3[ there exists Qo £ Qp(T) such that \Vq\ = iq , 
i.e. 

min < / H(x,\V\)dx : divU T, 1 _ / , q 

This shows that (j4.3j) holds true and that the infimum in the left-hand side is indeed a 
minimum. 

The relation between minimizers of the two problems is an easy consequence of the 
previous constructions. □ 

Remark 4.6. Similar Lagrangian formulations have been studied in connection with trans- 
port problems involving concave costs, like for example problems where to move a mass m 
of a length i costs m a I (0 < a < 1). For these the reader is referred to the monograph [3], 
as well as to the papers [29j EH] ■ 



Appendix A. Decompositions of acyclic flat 1-currents 

A.l. Definitions and links to vector fields. The classical references which we use for 
currents are [19\ 12 lj . We translate however all results in the language of the previous 
sections. In the following, by f2 C M. N we still denote the closure of an open bounded 
connected set, having smooth boundary. 
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Definition A.l. A 0-current on is a distribution on 0, in the usual sense. A 1- 
current on O is a vector valued distribution on 0,. The relevant duality is the one with 
1-forms uj{x) = Y2i=i dx l having smooth coefficients, i.e. C0i € C°°(Q). We denote by 
C°°(Q, A 1 R N ) the space of such forms. More generally, a fc-current is an element in the 
dual of smooth £;-forms C°°(ft, A k R N ). 

The above definition automatically gives the space of currents a natual weak topology, 
defined via the duality with smooth forms. For any current there is a natural definition of 
boundary. 

Definition A. 2. If J is a /c-current on f2, then we can define its boundary dl to be the 
[k — l)-current on Q which satisfies 

(dl,(p) = -(I, dip), for all p G C°°(n, A^V), 

where d is the exterior derivative. For example, if k = 1 we must take ip £ C°°(f2) and dip 
is the 1-form J2i d Xi <fi dxi . 

Definition A. 3. Let / be a /c-current, then its mass is 

M(J) = sup{ \{I,cj)\ : u e C°°(n,A k R N ), Bup||w(x)|| < lj , 

where the norm ||w|| for an alternating /c-tensor is defined as 

||w|| = sup{(w,e) : e unit simple k — vector}. 

For k = 1 this coincides with the usual norm \\cu\\ = W 'uf + . . . + uj^. 

We will be interested just in 1-currents which are distributions of order (i.e. vector 
valued Radon measures) and in their boundaries (which are scalar distributions of order 
1). For them, a comment is in order. 

Remark A. 4. Finite mass 1-currents can be identified with vector-valued Radon measures 
as follows. To every smooth 1-form io we may associate naturally a vector field X w := 
(wi, . . . , wjv). We can then write for a 1-current / 

J X u dl := (I,u)). 

Since C°° is dense in C°, the resulting linear functional on smooth vector fields can be 
identified via Hahn-Banach theorem to a unique linear functional on C° vector fields. The 
latter is indeed a vector-valued Radon measure by Riesz representation theorem. Since 
ll^ 2 -)!! = II^X^)!! by the above definition, we automatically obtain that 

M(J) = / d\I\, 
Jn 

i.e. the mass equals the total variation of I regarded as a Radon measure. The same 
reasoning can be applied to O-currents of finite mass, by identifying them with scalar 
Radon measure^. 

^See also "Distributions representable by integration" in |19l 4.1.7] 
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Definition A. 5 (variation of a current). Let A be a /c-current with k G {0, 1}. Then we 
may define the variation measure of A in the usual sense, by identifying A with a Radon 
measure as in Remark IA.4I Thus for a Borel set E we define 



(j, A (E) := sup < V / 

U=i M 



dA 



: Ei form a Borel partition of E 



An equivalent way of defining fiA would be as the infimum of all measures fi such that 
(A,oj) < J Q \\uj\\dfi for all smooth 1-forms u. 

We recall that a /c— current T is said to be normal ifQ 

M(T) + M(#T) < +oo. 

We now define flat currents, a class useful for its closure properties. 

Definition A. 6. We define the flat norm of a /c-current A as follows 

¥(A) = inf {M(,4 - dl) + M(J) : / is a (fc + l)-current with M(J) < oo} . 

Then the space of /?ai k-currents is defined as the completion of normal k— currents in the 
flat norm. 

Flat currents of finite mass have the following characterization, which will be exploited 
in the sequel. 

Lemma A. 7. Let T be a flat k— current of finite mass. Then T is a flat current if and 
only if there exists a sequence of normal k— currents {T n } n ^ such that 

lim M(T n - T) = 0. 

n— >oo 

Proof. This is a standard fact, we provide a proof for the ease of completeness. By definition 
of flat convergence, there exists a sequence {/ n } nS N of normal k— currents and a sequence 
{y n } ne N of (A; + 1)— currents such that 



lim 

n— >oo 



M(J -I n - 3Y n ) + M(Y n ) - - 

n 



< lim ¥(I - I n ) = 0. 

n— >oo 

We then set T n = I n + dY n , which by construction is a A;— current and of course 

lim M(T - T n ) = 0, 

thanks to the previous. Moreover, this is a normal current, since by triangular inequality 
we have 

M(T n ) < M(J„) + U(dY n ) < 2M(I n ) + M(T - /„ - dY n ) + M(T) < +oo, 
thanks to the fact that T has finite mass, and also 

U(8T n ) = M(dl n ) < +oo, 
since d(dY n ) = 0. This concludes the proof. 



For k = 0, we define dA — and thus the condition on dA can be omitted. 
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The converse implication is even simpler: by definition of flat norm, we directly have 

F(T - T n ) < M(T — T n ), 
which concludes the proof. □ 

A significant instance of flat 1— currents with finite mass is given by L 1 vector fields. We 
provide a small proof of this elementary fact. 

Lemma A. 8. Given V € L l (£l\M. N ), we naturally associate to it the 1 — current Iy of finite 
mass, defined by 

(Jy,w)=V f ujiVidx:= [ u(y), for every uj E C°°(0, A 1 !^). 
i=1 Jn Jn 

This has compact support contained in Q and M(7y) = ||V||^i. Moreover, Iy is a flat 
current. 

Proof. We just prove that Iy is a flat current, the first statement being straightforward. 
At this aim, we will use the characterization of Lemma IA.7I and we will construct the 
approximanting currents by convolution. For every e<l, we define 

Sl £ = {ieSl : dist(a?,dfi) > 2e}, 

then we take a standard convolution kernel g € C^°, supported on the ball {x : \x\ < 1} 
and we define 

N „ ( x \ _ ^ to,N 



Q e {x) = e Q\-)i xe 
We also set 

V e := (V-l^)*Qe, 

where 1e stands for the characteristic function of a set E. Consequently, we define I £ := Iy e 
and observe that V £ (and thus I e ) has compact support contained in £1. From the mass 
estimate and Holder inequality we obtain that masses are equi-bounded, since 

M(/ e ) < HKIIli < \\q\\l°°\\V\\ l i < C\\V\\ L i. 

The boundedness of dl £ follows a similar route: since V £ have compact support (strictly 
contained) in fi, we have 



1(01^)1 = 1(1^)1 = 

d\\V £ <p dx 



j 



V £ ■ Vip dx 

< ll div KllLi II^IIl 00 , 



so that setting C £ = |jdiv ||^i and passing to the supremum on tp, from the previous we 
obtain 

M(9J e ) <C £ < +oo. 

This implies that {I £ } £ >o is a sequence of normal currents. Moreover, the mass convergence 
M(/ £ — /) easily follows from the convergence of V £ to V in L 1 (il;IR Ar ). □ 
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Remark A. 9. With the previous notation, it is easily seen that the boundary dly corre- 
sponds to the distributional divergence of V, i.e. 

(dl v ,<p) = - [ V<^ • V dx for every ^ G C°°(0). 
Jn 

This distribution as well has compact support in f2. 

Definition A. 10. A 1-current I is called acyclic if, whenever we can write I = I\ + 12, 
with M(J) = M(/i) + M(J 2 ) and dh = 0, there must result h = 0. 

For I = Iy with ^ 6 L 1 , we have the correspondence with Definition 12.61 i- e - 
(A.l) V is acyclic Iy is acyclic. 

A. 2. Lipschitz curves as currents. We recall that we called £(f2) the space of equiv- 
alence classes of Lipschitz curves 7 : [0, 1] — > $7 (see the beginning of Section [4]), with the 
topology of uniform convergence. 

Here we remark that to each 7 G £(f2) we may associate a vector valued distribution, i.e. 
a 1-current, denoted by [7] and defined by requiring for all u G C°°(Q, A 1 M Ar ) 

p pi N pi 

<[ 7 W:= / w(7(*)) [V (*)] * = E / ^(7(*)K(*M- 

7 7 7o i=1 •/ 

Note that this expression is well-defined on C(Q), since the integral on the right is invariant 
under reparameterization. 

For 7 G £($7) we have M([ 7 ]) < ^(7) := Jq | 7 '(t)|(it with equality exactly when 7 has a 
representative which is injective for T^-almost every time. We are interested in a stronger 
requirement, namely that the curve does not even intersect itself, so that an injective 
representative exists (such curves are called "arcs" in |27|). We fix a notation for such 
classes of curves. 

Definition A. 11 (Arcs). We define £(£l) the subset of C(Cl) made of those classes of 
curves 7 : [0, 1] — > f2 which have an injective representative. 

Remark A. 12. One could think of the above-defined arcs as "acyclic curves", where a 
"cycle" can mean two things: 

• we can have a cycle in the parameterization, where a cycle would be represented by 
a curve satisfying (up to reparameterization) j(t) = 7(1 — t) and "inserting a cycle 
of type 7" in another curve 7 such that 7(^0) = 7(0) would result into the curve: 

f 7(2s) if a €[0, to/2] 

7(a) = { 7 (2s - t ) if s G [t /2, (t + l)/2] 
[7(2s-l) if a€[(t + l)/2,l]. 

• a curve which intersects itself (i.e. which has no injective parameterization) will 
instead give rise to a 1-current which is not acyclic, since it will have a reparame- 
terization containing an injectively parameterized loop. 
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The fact that in the decomposition of acyclic currents one restricts to using just arcs (for 
which neither type of cycle occurs) is then another natural consequence of the robustness 
of the acyclicity requirement. 

On 1-currents we consider the topology of distributions. The following result links the 
two topologies: 

Lemma A. 13. The map £(0) 3 7 [7] as defined above is continuous on the sublevels 
of the lenght functional £ : £(Q) — > M. 

Proof. Assuming that 7 j — > 7 and £{^i) < C we then obtain that 7 » converge uniformly. In 
particular they converge as distributions. □ 

A.3. Smirnov decomposition theorem. We are now in a position of stating the theorem 
on the decomposition of 1-currents due to Smirnov [36] and recenlty extended by Paolini 
and Stepanov in [27, 28J to metric spaces. The result which is relevant for us is contained 
in [36], though for its extension needed below we will use the techniques of |27| . 

Theorem A. 14 ( |27l 136] ). Suppose that I is a normal acyclic 1-current on fL Then there 
exists a finite Radon measure Q on the space of continuous curves on Q, concentrated on 
C(il), and such that the following decompositions of I are valid in the sense of distributions: 



(A.2) /= / [ 7 ]dQ( 7 ) and m = / A* [7 ]dQ( 7 ); 

Jc(ci) Jc(n) 

and 

(A.3) dl = d[y]dQ{i) and p dI = / ^1 dQ{i). 

Jc(n) Jc(n) 

We now note down some reformulations of the items present in the above theorem, in 
terms of measures and vector fields: 

• the total variation is the mass norm, i.e. /x/(0) = M(J) and /j,qi{Q) = M(5J); 

• if V is a L\ oc vector field, then ly has variation measure \ij v = \V\ ■ Jz? N ; 

• if p = p + — p- is the decomposition of a signed Radon measure into positive and 
negative part, then p p = p + + p_; 

• in particular for 7 € £(0) we have (AgM = ^ 7 (i) + ^7(0)- Since this measure has 
total variation 2 for all 7 , we can quantify the total mass of the above Q by means 
of the mass norm of the boundary of /. Namely, we have 

Q(c(n)) = l m (n) = ^p-, 

• for 7 6 C(fl), there holds //u = J^ 1 Llm( 7 ) is the arclenght measure of 7 ; 

• by expanding the definitions and comparing to Section [H we see that for I = Iy 
with V € L 1 (0) and for Q as in Theorem IA.141 there holds 

\v\ ■ sf N = m v = % Q ■ se N , v ■ j? n = 1 = i Q ■ x N 
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and 

div V = dly = (eo - ei)#Q. 

All these reformulations allow to translate Theorem IA. 141 in the case of I = Jy, with 
V € L p (f2), for p > 1. This is the content of the next result. Observe that since for a 
bounded Q we have the continuous inclusion L p (£l) L 1 (fi), we can restrict the result to 
p = 1 without loss of generality. 

Corollary A. 15 (Reformulation of Theorem I A. 14]) . Suppose that V € L 1 (f2) is an acyclic 
vector field such that div V is a signed Radon measure. Then there exists a finite Radon 
measure Q on the space of continuous curves on Q, concentrated on £(0), and such that 
the following decompositions of V are valid in the sense of distributions: 

(A.4) iq = V and i Q = \V\, 

and 

(A.5) - divV = (ei - e ) # Q and (div V)+ + (div V)_ = {e 1 + e )#Q. 

We are now lead to consider the extension of Theorem IA. 141 to the case where div V 
is not a Radon measure. This is necessary since in our problem divV E W~ 1,P {£1), so 
in general it is not a measure. Some non trivial issues arise in this case: first of all, the 
measure Q which decomposes I may not be finite, in general. 

Example A. 16. For p < j^zj, we consider an infinite sequence of small dipoles {(ai, bi)} i 
such that 

oo 

y~] |o, — 6 i | 7V ~P( iV ~ 1 ) < +oo, and Ai;A are disjoint, 

-i=i 

where as in Example 12.41 we used the notation 

D ai , bi = U Xl ,x')eR N : \ x '\ + \ Xl \<h^i\. 

If we consider the vector fields V ai & 4 as in Example 12.41 then the new vector field defined 
by V = Yh=i V ai,bi verifies 

P oo 

<C^^|a l -6,r-^ Ar - 1 )<+oo, 
LP i=l 

which implies T := - $>J G W ,_1 'P(n). By observing that oo = M(T) = J £(Q) dQ 

for any decomposing measure, we see that no finite measure Q can be found. On the other 
hand, a a— finite measure Q can be found, since each V ai ^ t can be separately decomposed 
with a measure Qi of mass 2 and the have disjoint supports. 

Example A. 17. The following is another version of Example IA.161 which exploits the 
Sobolev embedding theorem. Let us take q > N, then W 1,q (ri) can be identified with a 



w 



IP 

I Li-' 



E 

i=l 



V, 



ai,bj 
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space of functions which are Holder continuous of exponent a = 1 — N/q. Then we pick 
the following two curves (here £< 1) 



71 (t) 



1 



(cost, sin t) 



t 2 / a 

We define the distribution 

{T,<p) = 



and 



72 (t) 



1 



t 2 / a 



(cost, sini), t > 1. 



V(7i(*))-V(72(*)) ^GC7°°(0), 



then again this is an element of W 1,P (Q), since by Sobolev embedding [(/?]c°> a < Cfi IMIvy 1 . 1 ?' 
so that 



|<2»|< / \^m-vMt))\dt<cy\\ wl , 9(a) l \ 7l (t)- 72 (t)\ a dt 

= Ce a ||y||vvi,9(n), 



We then introduce the measure on paths Qt, defined by 



Q 



T 



5 



71 W72W 



(it. 



where for every t > 1, by 71 (i) 72 (t) we indicate the straight segment going from 71 (t) to 
72 (t). Observe that for every ip we have 



£(H) 



K7(0))-^(7(l))]dQr(7) 



^(7l(*))-^72(t)) dt={T,<p), 



and 



£(n) 



%)dQr(7) 



l7i(t)-72(*)|cft 



t 2 A 



eft 



a 



a 



00 a 
= e 

1 2-q 



< 00, 



while Qt is not finite, but just c-finite. 



The previous examples clarifies that we cannot hope to give a distributional meaning to 
the positive and negative parts of the divergence of V, and the good definition of div V as 
a distribution relies in general on some sort of "almost-cancellation" . Therefore the last 
no-cancellation requirement (div V)+ + (div V)_ = +en)-#Q of Smirnov's Theorem lA.141 
must be relaxed when we extend it to a larger class of Vs. 

We introduce a notion needed in the next theorem: 

Definition A. 18. Let A,B,C be 1-currents. We say that C is a subcurrent of A if 

M(A - C) + M(C) = M(A). 
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In this case we write C < A. If also M(B - C) + M(C) = M(B) we say that C is a common 
subcurrent of A and -B. We will use the notation 

SC(A,B) = {C is a 1- current : C < ^ and C < B}, 

to denote the set of common subcurrents. 

Let C G SC(A, B), then we say that C is a largest common subcurrent of A and -B if 

SC(A-C,B -C) = {0}, 

where by we denote the trivial 1— current such that (0, 99) = for every 1— form (p. In 
this case, we will use the notation C = A A B. 

Lemma A. 19. Let A,B be normal 1-currents. Then they admit a largest common sub- 
current C. 

Proof. The proof is exactly as in [27, Proposition 3.8] so we just indicate the main modifi- 
cations needed in our case. For every pair of currents (T, S), we set 

£(S,T) := sup{M(C) : C G SC(S,T)}. 

Let us define {Ay, B u , C u } ue n recursively as follows 

A\= A — Co, B\ = B — Co, 

C G SC(A , B ) such that M(C ) > , 
A„+\ = A v — C, B v+ \ = B v — C v 

where C v G SC{A V , B v ) such that M(C„) > Bv) for u>l. 

Then as in [27] we have that M(C) < £(A,^)/2 for all C G S'C^+i, B v+l ) and in 
particular by taking the supremum we can infer 

f( ^ +l) <«^<...<«^^0. 

From this we obtain that {AjveN and {By}y £ ^ are converging in mass norm and by setting 

A' : = lim A v and 5' := lim B v , 

we obtain A' < A v and B' < B v , for every v G N. This implies that any C G SC(A', B') 
also belongs to SC(A U ,B V ) and by definition we have M(C') < £(ylj,, — > 0, thus we 
get C = 0, i.e. the trivial current is the only common subcurrent of A' and B' . We then 
observe that A — A' = B — B', since A v _\ — A u = B v -\ — B u for v > 1 by construction, so 
we have 

A — A' = lim L4 - A u ] = lim [C +A\- A v \ 

= lim [C + {A 1 -A 2 ) + --- + - A/)] 

= lim [C + {B x - B 2 ) + ■ ■ ■ + (By_! - By)} 

V — ^00 

= lim [B — By] = B - B'. 
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We can finally set C = A - A' = B - B' € SC(A, B) and obtain SC(A -C,B - C) = 
SC{A' , B') = {0}, i.e. C is a largest common subcurrent of A and B. □ 

Theorem A. 20 (Smirnov Theorem for flat acyclic currents). Let I be an acyclic flat 1- 
current, having finite mass. Then I is decomposable in curves, i.e. there exists a a- finite 
positive measure Q on curves such that Q is concentrated on £(£l) and 



(A.6) /= / h]dQ( 7 ). 

Jc(n) 

Moreover, we have 

(A.7) m = / dQin) 

Jc(n) 

and 

(A.8) 81= [ 9[ 7 ]dg( 7 ). 

Jc(n) 

Proof. We start observing that thanks to Lemma IA.71 our current / is approximable by 
normal currents {Jj}j g N in the mass norm, i.e. limM(J — Jj) = and M( Jj) +M(<9Jj) < oo. 
We construct a first decomposition of / using this information. 

Step 1. (I can be decomposed in normal currents). We may assume that {JijieN is a 
Cauchy sequence in the mass norm, i.e. that M( Jj_i — Jj) < 2~ l for all i. We then define 

-^o = Jioi 

h = Ji +k — Jio+k-i for k > 1. 

It then follows that I = Y,kh and J2k M (^k) < 2M(J) if i £ N is chosen large enough. 
We then modify the Ik into a better decomposition Ik as follows: 

4+1 = 4+1 - 4+1 A 4 for k > 0. 

where we recall that 4+1 A 4 denotes a largest common subcurrent of 4+1 and 4) which 
is well-defined thanks to Lemma IA.19I From these definitions it follows that we have the 
decomposition 

7 = ^4, M(/) = ^M(4)- 

fc>0 k>0 

Moreover the currents 4 are still normal since each of them is a finite sum of subcurrents 
of normal currents. 

Step 2. (Use of Smirnov's result for normal currents). By applying Theorem I A. 141 to 
each of the 4 we obtain a sequence of measures Qk on the set C(Q) concentrated on £(£1), 
such that Qk decomposes 4 in the sense of Theorem IA.14I From the previous step we 
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obtain then a decomposition for / as follows: 

(A.9) / = E j * = E/ bWkb), 

k>0 k>0^ C ^ 

and 

(A.10) M(J) = ^M(/ fc ) = E / - 

Step 3. (Construction of a a-finite measure). We may then define a measure on £(Q) by 
setting 

Q(A) = ^ Qfc(-4), for every Borel set A C 

fc>0 

where the series on the right has to be intended as ^fc>o Qk{A) = lim^oo Ylk=o Qk{A) and 
this limit exists by monotonicity (it is in fact a supremum), since the measures {Qk}ke~N 
are positive. At first, we verify that Q is indeed a measure: we have 

oo 

Q(0) = Qk{®) = 0, as Q fc (0) = 0, for every k G N, 

fc=0 

and also for any collection {^i}^ of disjoint Borel sets in £(£l) we get 

(oo \ V / oo \ oooo 

i=l / fc=0 Vi=l / fc=0 8=1 

oo oo oo 

= EE^(^) = E^)> 

i=l fc=0 i=0 

where we have the right to exchange the two series, since each term Qk(Ai) is positive. 
This finally shows that Q is a measure on C(Q). We also observe that 

oo 

Q (£(fi) \ £(fi)) = Qk \ = 0, 

fc=0 

as each Qk is concentrated on £(f2). 

In order to verifiy that Q is a— finite, for all h G Z we define the Borel sets -B^ := {7 G 
: £(7) G [2 h , 2 h+1 [}, which form a partition of Observe that the fact that B h is 

a Borel set easily follows from the lower semicontinuity of the length functional. Moreover, 
for fixed h G Z there holds 

Q(^) = ^Q fc (B ft ) < 2- h ^M(I fe ) = 2- ft M(/) < +00. 

fc>0 fe>0 

This finally shows that Q is a— finite. The fact that Q verifies (1A.6|) and (IA.7jl follows 
directly from (TOj) and (lAlfl| . while (EOjl is a consequence of (|A.6j) , which suffices to test 
on special forms to = d<\>. □ 
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Remark A.21. The proof of the Theorem IA.20I extends to general currents on metric 
spaces, following [27] (see also [28] and [36]). 
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